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Abstract 

Using density functional theory combined with quasi-harmonic approximation, we investigate 
the thermal expansion behaviors of three different types (ct, /3, and 7 ) of graphyne which is a 
two-dimensional carbon allotrope composed of sp and sp^ bonds. For each type of graphyne, we 
obtain the temperature dependent area variation by minimizing its free energy as a function of 
temperature, which is calculated by considering all the phonon modes in the whole Brillouin zone. 
We find that all three types of graphyne exhibit negative in-plane thermal expansion up to T ^ 
1000 K. The observed in-plane thermal contraction can be attributed partially to the ripple effect, 
similarly in graphene. The ripple effect itself, however, is not sufficient to explain anomalously 
larger thermal contraction found in graphyne than in graphene. Our deliberate analysis on the 
phonon modes observed in graphyne enables us to discover another source causing such thermal 
expansion anomaly. We find that there are particular phonon modes with frequencies around a 
few hundreds of cm“^ existing exclusively in graphyne that may fill empty spaces resulting in area 
reduction. These modes are identified as “rigid unit modes” corresponding to the libration of each 
rigid unit composed of sp^ bonds. 
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I. INTRODUCTION 


Various carbon allotropes based on sp^ hybridization, snch as fnllerenes, [Ij nanotnbes, [2] 
and graphene, [3] have been intensively studied during last a few decades, since their first 
discoveries. Their unique properties are basically related to the hexagonal geometry and the 
Dirac cones in the electronic structure of graphene, [3Hi which is a base material to form 
other sp^-bonded carbon allotropes. To make nse of graphene in future nanoelectronics, 
however, it is indispensable to open an energy gap at the Fermi level E^. Althongh lots 
of attempts have been made to open a band gap, [M2] nothing has yet come along for 
practical applications. As a separate approach, sundry efforts have been made to search 
for other 2D allotropes intrinsically possessing an energy gap at Ep to be used as graphene 
substitutes in electronic applications. [13] 

A new 2D layered carbon allotrope called graphyne was proposed in 1987. [H] It can be 
geometrically generated from graphene by inserting a single acetylenic linkage (AL) consist¬ 
ing of two sp-hybridized carbon atoms into a bond in graphene composed of sp^-hybridized 
carbon atoms only. Since there are a number of ways of inserting ALs, many different 
types of graphyne can be formed. Among them were mainly considered three types with 
a hexagonal symmetry, classified by the number of ALs attached to each sp^-hybridized or 
triply-coordinated carbon atom: «-, /3-, and y-graphyne. Fig. [^a-c) shows their equilib- 
rinm structnres. In a-graphyne shown in Fig. [^a), there are three ALs connected to each 
sp^-hybridized carbon atom, whereas in /5-graphyne (y-graphyne), every sp^-bonded atom 
has only two ALs (one AL) as shown in Fig. [^b) and (c), respectively. 

Theoretical studies reported that both a- and /3-graphyne structnres possess Dirac cones 
near Ep in its Brillouin zone (BZ) representing their semi-metallic characteristics like in 
graphene. [ISH2I] It was also reported that y-graphyne, which is by far the most stable 
graphyne, is interestingly an intrinsic semiconductor with a direct band gap of 0.46 eV at M 
points. The existence of the band gap was ascribed by the Peierls instability resulting in the 
wave fnnction localization at the triple bonds. [I8] It is, therefore, expected that y-graphyne 
could be made use of as a semiconducting component in future nanoelectronics. 

Although various research groups have endeavored to synthesize them, graphyne has not 
yet been realized in a crystalline form. Nevertheless, it was already reported that graphyne 
sub units that resemble y-graphyne were successfully synthesized. [22H21] Moreover, one of 
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FIG. 1. Equilibrium structures of (a) (a-graphyne, (b) /3-graphyne and (c) y-graphyne. Respective 
hexagonal unit cells are depicted by dashed parallelograms. Nonequivalent bonds in each type of 
graphyne are denoted by di, ^2, and df. di and ^2, which are single-bond like, represent a bond 
between an sp^-bonded (or triply-coordinated) and an sp-bonded (or doubly-coordinated) atoms, 
and one between two triply-coordinated atoms, respectively, whereas dt designates a triple bond 
between two sp atoms in every acetylenic linkage, (d) Cohesive energy comparison of three types 
of graphyne and graphdiyne with respect to graphene as a function of Nsp/N^ where Nsp and N 
are the number of sp-bonded or doubly-coordinated atoms and N and the total number of atoms 
in the unit cell, respectively. 
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their cousin structures with two ALs in a sequence instead of a single AL in 7 -graphyne, 
called graphdiyne, was also reported to be discovered in the form of sheet, [2S] as well as in a 
tube form. [SH] According to our calculation as well as other studies, [TTfEH] /3-graphyne and 
7 -graphyne are even more stable than graphdiyne, as shown in Fig. Bd). and thus expected 
to be synthesized in near future. 

Compared to the structural and electronic properties of a-, (5-, and 7 -graphyne that 
have been intensively studied as discussed above, there have been not many studies on 
their thermal and thermal expansion (TE) behaviors, [201 EH EO] which are also necessary 
properties to be disclosed in advance for them to be utilized as various devices in future 
nanoelectronics. 

It is well known that most materials expand thermally with a very different material- 
dependent linear TE coefficient (TEC) of 10“^ ~ 10“^ mainly due to the atomic vibration 
under asymmetric potential. In almost all devices and equipment, such a TE usually be¬ 
comes a critical problem deteriorating their performance and lifetime. There have been, 
on the other hand, several reports revealing materials with a negative TEC (NTEC), such 
as graphite, EH carbon nanotubes, [HUES! zirconium tungstate (ZrW 208 ), [31| perovskite 
oxides, [HSHSZl and Invar alloys, [SHlEni originating due to a ripple effect, bending/twisting 
modes, rigid unit modes (RUMs), jlHlIlI] atomic radius contraction, [12] or a magnetovolume 
effects. [l3| Researchers have made various attempts to make composites containing NTE 
materials, such as ZrW 208 [lH - 116] and /3-eucryptite m , as a TE compensator, expecting 
them to be high-performance composites with near zero TE. It was indeed shown that TE 
can be suppressed in molten aluminum alloy with NTE manganese antiperovskites. [48] 

In this letter, we report a first discovery of RUMs in two-dimensional graphitic carbon 
systems, which are responsible for a TE anomaly of NTE in sp-sp^ hybridized carbon sys¬ 
tems, similarly observed in 2D perovskites. H The thermal properties and TE behaviors 
of a-, (5-, and 7 -graphyne were investigated by first-principles calculations of temperature- 
and volume-dependent free energies from the phonon properties. Similar to graphene, all of 
three configurations contract with temperature up to T~1000 K, but the absolute values of 
their TECs are much larger than that of graphene especially at low temperatures. Based on 
our thorough investigation on their phonon modes, we attributed such thermal contraction 
behaviors to not only the ripple effect as in graphene, but also RUMs corresponding to 
libration motions of rigid units composed of sp^-bonds only. 
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II. COMPUTATIONAL DETAILS 


To investigate the TE behaviors of three types of graphyne, we first obtained their equi¬ 
librium structures by carrying out first-principles calculations within the density functional 
theory (DFT), [50] as implemented in Vienna ab initio simulation package (VAST), [5T],|52| 
while adjusting three different carbon-carbon bonds, and thus the lattice constant of each 


graphyne as discussed in Sect. Ill, Projector augmented wave potentials [sniEii was em¬ 
ployed to describe the valence electrons, and the electronic wave functions were expanded by 
a plane wave basis set with the cutoff energy of 450 eV. The exchange-correlation functional 
was treated with the Perdew-Burke-Ernzerhof (PBE) parameterization (53] of the gener¬ 
alized gradient approximation (GGA). To mimic single-layered graphyne, we introduced a 
vacuum region with 15 A along the c-axis perpendicular to the sheet. The BZ was sampled 
using a P-centered 10x10x1 fc-grid for the primitive unit cell of each type of graphyne. Al¬ 
though the primitive cell size of /3-graphyne is larger than the other types, the same fc-grid 
was employed for precise calculation for /3-graphyne because there were earlier studies on 
/3-graphyne showing inconsistent results on structural stability. [201 E] 

The phonon dispersion relations were computed by applying the finite displacement 
method (ED) [56l [57] for 3x3x1 supercells. The corresponding reduced BZs were sampled 
by 3x3x1 fc-point meshes. In each supercell, we solved the secular equation of dynamical 
matrix constructed at every wave vector q from the force constant matrices computed under 
ED to obtain the phonon dispersion relation. To investigate thermal properties, we em¬ 
ployed quasi-harmonic approximation (QHA), [58l [5^ in which the Helmholtz free energy 
of a 2D graphyne sheet was calculated by 


F{T, A) = U{ A) + - 


^ E'^1'“ + EIn |l - exp } 

q,n q,n ^ ^ ^ 


( 1 ) 


where h, and ks are the reduced Planck constant, and the Boltzmann constant, respectively; 
and is the phonon frequency with the wave vectors q and the mode indices n. In Eq. Q, 
the first term U(A) is the system internal energy with a constant area A at T = 0; the second 
term represents the vibrational zero-point energy of the lattice; and the last term corresponds 
to the phonon contribution to the free energy in QHA. [58l[59] U (A) and a;q^„ were evaluated 
at 12 different area points. F{T,A) was fitted to the integral form of the Vinet equation 
of state (HOI to obtain the minimum values of the thermodynamic functions with respect to 
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the area, and thus the equilibrium area, A{T), as a function of temperature. TE behaviors 
of three types of graphyne sheets were explored by calculating the area TEC as{T) as a 
function of T, 


as{T) 


1 dA{T) 
A{T) dT 


( 2 ) 


III. RESULTS AND DISCUSSION 

We first searched for the equilibrium structures of three kinds of graphyne, which are 
indispensable to their phonon dispersion relations and thermal properties. Graphyne can 
be made by replacing all or some bonds in graphene by an AL (-C=C-). Among numerous 
different types of graphyne produced in this way, we only focused on three types of gra¬ 
phyne, a, /3, and 7 , which are highly symmetric in a hexagonal lattice relative to the other 
types. The a-graphyne can be made by a complete substitution for graphitic bonds, while 
the /3-graphyne and 7 -graphyne can be formed by two-thirds and one-third substitutions, re¬ 
spectively. In these types of graphyne, there are at most three nonequivalent carbon-carbon 
bonds, di, ^ 2 , and dj. Single-bond like di and d 2 represent bonds from a sp^ atom to a 
neighboring sp and sp^ atom, respectively, while a triple-bond dt describes a bond between 
two neighboring sp atoms forming an AL together with two di bonds. 

Eor each graphyne, we completely scanned the whole energy surface while adjusting the 
bond length of three nonequivalent carbon-carbon bonds, di, d 2 , and dt keeping its hexagonal 
symmetry. Note that there are no d 2 bonds in «-graphyne, and the lattice constant of 
each unit cell is uniquely determined for each set of given di, d 2 , and dj. The equilibrium 
structures of three types of graphyne were determined at the minimum energy point in 
each energy surface, and shown in Eig. [^a-c), where their hexagonal lattices were depicted 
by dashed parallelograms representing their unit cells. Two nonequivalent bonds in the 
«-graphyne were calculated to be di = 1.397 A and dt = 1.230 A with lattice constant of 
Ueq = 6.969 A. Eor the /5-graphyne ( 7 -graphyne), on the other hand, its calculated lattice 
constant and three nonequivalent bond lengths are Ueq = 9.480 A (6.890 A), di = 1.389 A 
(1.408 A), d2 = 1.458 A (1.426 A), and dt = 1.232 A (1.223 A). Our determined structural 
parameters are in good agreement with other studies. [HI | 2 T] 

To compare the structural stability of three types of graphyne, we evaluated their cohesive 
energies with respect to that of graphene as a function of Ngp/N, where Ngp and N are the 
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FIG. 2. The phonon dispersion relations of (a) a-graphyne, (b) /3-graphyne, and (c) y-graphyne 
along the high symmetry lines in 2D hexagonal Brillouin zones. 



number of sp-bonded or doubly-coordinated atoms, and N and the total number of atoms 
in the unit cell, respectively. As displayed in Fig. [^d). these types of graphyne are only 
0.6 —0.8 eV/atom less stable than graphene. This tendency agrees well with quantum Monte 
Carlo calculations [2M9|. Especially, the y-graphyne and /3-graphyne are even more stable 
than graphdiyne, which was already reported to be synthesized, 123 Ei expecting to be 
synthesized as well. 

Using 3 X 3 X 1-repeated supercells of the equilibrium structures, we calculated the phonon 
dispersion relations of all three types of graphyne as shown in Fig. and identified all the 
phonon branches by analyzing their corresponding eigenvectors. We found in all cases that 
there are three acoustic phonon branches below u ~ 400 cm“\ two in-plane modes with a 
linear dispersion and one quadratic out-of-plane mode near F point. We did not find any 
imaginary frequencies in all the types of graphyne indicating their structural stability. Note 
that there were inconsistent results for /5-graphyne in previous studies showing imaginary 
phonon frequencies iza or phonon mode anomalies, [ 20 ] which were not observed in our 
calculation. We believe that such inconsistency was due to use of an incompletely relaxed 
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structure in a small supercell size. On the other hand, the calculated phonon dispersion 
relations of cr-graphyne and 7 -graphyne are in good agreement with other studies. [201 El] 

Compared to graphene with its maximum frequency of .^1600 cm“^, | 6 T] graphyne ex¬ 
hibits even higher phonon branches dispersed at ^ 2000 cm“^. These high frequency 
modes correspond mainly to the stretching vibration of the triple bonds dt with larger force 
constants than single-bond-like di and ^ 2 , which generate graphene-like phonon modes at 
u ^ 1500 cm“^ below a large frequency gap of Au ~ 500 cm“\ as displayed in Fig. 
The existence of such a large frequency gap implies that the vibration of dt bond stretching 
modes is strongly protected from being scattered by the graphene-like modes by the energy 
conservation law. [62| Another phonon gap of Au 32, 78, or 67 cm“^, was also observed 
near 500, 750, or 1000 cm“^ for a-, (3-, or 7 -graphyne, respectively. We found that the 
di and ^2 stretching modes are dispersed above the gap, while all the bending modes includ¬ 
ing optical flexural modes are settled below. Note there is an exception in the a-graphyne 
as displayed in Fig. [^a) that there are two flat modes identified as optical flexural modes 
located at zz ~ 602 and 622 cm“^ above the gap. Therefore, in the a-graphyne, these optical 
flexural phonons can be scattered by the stretching phonons, while in the other types of 
graphyne, they are located at similar frequencies, but below their frequency gap, as shown 
in Fig. [^b) and (c), which protects them from being scattered by their corresponding di 
and 6^2 stretching modes. 

To evaluate the Helmholtz free energy, we used Eq. Q with the phonon dispersion rela¬ 
tions cuq^n calculated at 12 different area points, based on QHA. For each type of graphyne, 
the calculated Helmholtz free energy was plotted as a function of temperature and area 
in a color-coded map as displayed in Fig. Note that color-coded energy values clearly 
represent that the 7 -graphyne is more stable than the other two types. The temperature 
dependence of the equilibrium area A(T) was obtained by minimizing the thermodynamic 
functions with respect to the area as denoted with the black solid line on each color-coded 
map in Fig. 

We explored the TE behaviors of three types of graphyne as well as graphene for com¬ 
parison by evaluating the temperature dependence of change in each system’s area with 
respect to the zero-temperature area, AA/Aq = [A{T) — Ao]/Ao, and its area TEC as{T) 
given by Eq. ([^. It was found that similar to graphene, all three types of graphyne contract 
as temperature increases up to T ~ 1000 K, as shown in Fig. |^a). Such an area shrink 
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FIG. 3. (Color online) Color-coded maps of Helmholtz free energy of (a) cr-graphyne, (b) (3- 
graphyne, and (c) y-graphyne as a function of temperature and area calculated using Eq. The 
black solid line on each plot indicates minimum value of free energy at each given temperature 
with respect to the area, and thus the equilibrium area for a given temperature. The color-coded 
free energy values are given in the color bar below. 

observed in a 2D sheet was attributed to its flexural modes, which are mainly responsible 
for the ripple effect, [63] as in graphene, [SOI ED and in graphyne. [SUED] We confirmed the 
contribution of flexural modes of graphyne to NTE by evaluating their corresponding mode 
Griineisen parameters, 7 q_„, defined as 

_ dlog _ V dujq^n 

aiogG “ dV ' 

to be negative. 

Using Eq. (§ , we also calculated the TECs of three types of graphyne as well as graphene. 
Compared to our calculated TEC values of graphene, which are in good agreement with 
previous studies, j30l |6T] it was observed that graphyne exhibits significantly lower TEC 
values than graphene as shown in Fig. ib). which cannot be fully understood only by 
the flexural modes responsible for the ripple effect. To scrutinize such large differences, 
we thoroughly inspected the real-space vibrations of all the phonon modes by visualizing 
their corresponding eigenvectors together with mode Griineisen parameters at different wave 
vectors q and mode indices n, and found intriguing unexpected vibrational modes, which 
we classified as RUMs. BS| It has already been known that the RUMs and quasi RUMs 
(qRUMs) involving a little distortion of rigid polyhedra are major causes of NTE occurred 
in an oxide framework consisting of rigid polyhedra, such as MO4 and MOe, where M and 
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FIG. 4. (a) Temperature dependence of the ratio of area expansion to the reference area at T = 0 K, 
AA/Aq in the unit of percentage %, and (b) the corresponding area thermal expansion coefficient, 
as{T)^ defined by Eq. in the unit of 10“^ of a-graphyne (dotted line), /3-graphyne (dashed 
line), and y-graphyne (solid line) as well as of graphene (dash-dotted line) for comparison. 

O are a metal cation and oxygen, [64l [65] although it was reported that there is no simple 
correlation between NTE and RUMs. [ 66 ] 

Fig. a a-c) shows real-space visualization of typical RUMs (in the right) observed at 
u 382, 277, and 249 cm“^ in a-, jS-, and 7 -graphyne, respectively, in tandem with their 
corresponding equilibrium configurations (in the left). The corresponding rigid units are 
denoted respectively by gray polygons of a triangle, a rectangle, and a hexagon, each of 
which encloses the sp^ carbon atoms connected with di and ^2 bonds. During the libration 
mode of the rigid unit in each graphyne, every AL becomes bended and thus the dt bonds 
are lengthened as shown in Fig. [^a-c). We found that such lengthened bonds can be 
compensated by filling an empty space or reducing the area size. To clarify the correlation 
between these RUMs and NTE in the graphyne systems and to analyze the size reduction 
quantitatively, we calculated the relative energy of each type of graphyne as a function of the 
area change AA/Aq with respect to the zero-temperature area Aq, and the rotation angle 
(p of each rigid unit. Fig. |^d-f) displays the color-coded energy map, in which the area 
reduction was estimated by searching for the minimum values of the relative energy AE 
for a given cp, as depicted with cyan-colored line with dots. Since the maximum rotation 
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FIG. 5. Real-space visualization of typical rigid unit modes observed in (a) cr-graphyne (b) /3- 
graphyne, and (c) y-graphyne respectively with v 382, 277 and 249 cm“^. Each local rigid 
unit is represented by gray polygon (triangle, square, or hexagon in (a), (b), or (c), respectively) 
surrounding carbon atoms with d\ and d 2 bonds. Each equilibrium structure is depicted in the 
left box, while the right box shows a snapshot of rigid units rotated counterclockwise corresponding 
to the selected RUM. Relative energy map color-coded as a function of AA/Aq and rotation angle p 
of each rigid unit calculated for (d) (a-graphyne, (e) /3-graphyne, and (f) y-graphyne. Cyan-colored 
solid line with dots indicates the area contraction determined by the minimum energy for given p. 
Energy values are given in terms of the equilibrium energy of each graphyne. 
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angle may increase with temperature, RUMs as well as the ripple effects may contribute 
to the area reduction. The mode Griineisen parameters corresponding the RUMs was also 
calculated to be negative as similarly in the flexural modes contributing to the ripple effects. 


IV. CONCLUSIONS 

We presented the phonon dispersion relations and the thermal expansion behaviors of a-, 
jS- and 7 -graphyne using the density functional theory. Our calculated phonon dispersion 
relations showed that there were no imaginary frequencies implying their structural stability. 
Their high-frequency modes at ~ 2000 cm“^ correspond to the stretching modes of the 
triple bonds dt, which are protected by the large frequency gap from being scattered by 
graphene-like modes of di and d 2 bonds. Similarly, the stretching modes of di and d 2 bonds 
can also be protected from being scattered by their bending modes, since the latter modes are 
separately located below another frequency gap, with an exception in a-graphyne. Our quasi 
harmonic approximation calculations revealed that all types of graphyne exhibit negative 
thermal expansion as in other 2D materials, but they have much larger values in negative 
thermal expansion coefficients than graphene. This significant discrepancy was resolved 
by identifying the rigid unit modes, which were firstly observed in 2D carbon graphitic 
structures. We found that not only the ripple effect originating from the flexural modes, 
but also the rigid unit modes are strongly responsible for the negative thermal expansion in 
graphyne. 
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